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Title:  Time-resolved  Near  Infrared  Optical  Tomography 
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CUNY  RF#  447349 


PI:  R.R.  Alfano. 

Researchers:  A.Ya.Polishchuk(Theory),  F.  Liu(Expt.),  B.B.Das(Expt.),  J.Dolne  (exp) 
Grad.  Students:  M.  Zevallos,  Y.X.  Shan. 

Goals: 

To  develop  theoretical  models  and  experimental  procedures  based  on  the 
detection  of  early  arriving  light  to  characterize  inhomogeneities  in  turbid  media  in 
multiple  scattering  regime  with  spatial  resolutions  on  the  order  of  transport  mean  free 
paths  (1-2  mm  for  biological  tissues  and  highly  scattering  media). 

Objectives: 

During  the  reported  period  the  following  objectives  were  pursued: 

1.  to  investigate  the  limits  of  applicability  of  the  diffusion  approximation  (DA) 
for  source  detector  distances  large  compared  with  photon  transport  mean  free  paths 
(where  the  DA  is  commonly  believed  to  hold),  and  show  that  for  higher  frequency  modes 
of  photon  migration  and/or  when  absorption  is  about  -10%  of  the  extinction,  the  DA 
becomes  a  progressively  approach  as  the  source  detector  distance  increases. 

2.  to  develop  a  new  approach  free  from  this  drawbacks  of  DA. 

3.  to  clarify  the  use  of  non-Euclidean  diffusion  (NED)  equation  especially 
tailored  to  describe  non-diffusive  features  of  photon  migration  in  highly  forward 
scattering  media  and  to  show  the  relationship  between  NED  and  the  more  general  theory 
of  radiative  transfer  (Boltzman  equation). 

4.  to  develop  a  path  integral  solution  of  the  NED  equation  in  the  case  of 
macroscopically  inhomogeneous  turbid  media. 

5.  to  perform  a  series  of  time-resolved  experiments  to  check  on  the  predictions  of 
the  NED  equation  against  the  intensity  temporal  profiles  of  scattered  light  measured  at 
various  distances  between  the  source  and  the  detector. 

6.  to  investigate  the  non-diffusive  nature  of  photon  paths  in  turbid  media  by 
analyzing  the  experimental  results  using  the  NED  equation  and  find  out  the  potential  of 
this  equation  to  improve  the  resolution  in  future  optical  tomographic  reconstruction. 
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Approach: 

A  close  interaction  between  theory  and  experiment  underlies  our  efforts  to 
understand  fundamental  features  of  photon  migration  in  turbid  media,  in  order  to  improve 
the  resolution  in  optical  tomography. 

The  experimental  part  of  the  program  was  based  on  the  measurement  of  temporal 
intensity  profiles  of  the  scattered  light  in  turbid  media  using  100  fs  laser  pulses  and  a 
streak  camera  with  10  ps  resolution.  The  turbid  media  used  were  biological  samples  like 
chicken  breast  tissue  as  well  as  model  random  media  with  monitored  scattering  and 
absorption  properties:  aqueous  intralipid  solutions  and  polystyrene  bead  suspensions. 

In  the  theoretical  part  of  the  program  numerical  solution  of  the  radiative  transfer 
equation  was  used  to  investigate  various  generalizations  of  the  diffusion  equation 
including  a  new  scalar-wave  diffusion  equation(SWDE)  introduced  in  the  course  of  this 
study.  This  part  of  the  program  was  aimed  at  achieving  a  better  description  of  photon 
migration  in  isotropically  scattering  media  at  large  distances  for  high  frequency  diffusive 
modes  and/or  with  absorption  exceeding  -10%  of  the  total  extinction. 

The  concept  of  photon  diffusion  on  the  velocity  sphere  and  non-Euclidean  path 
integral  technique  both  introduced  in  the  frame  of  this  grant  underlie  our  analysis  of non- 
diffusive  features  of  photon  migration  in  multiple  scattering  regime  in  highly  forward 
scattering  media.  This  approach  was  used  to  interpret  time  resolved  experimental  data, 
which  demonstrate  the  existence  of  ‘preferred’  photon  paths  on  scales  of  the  order  1,  (or 
smaller)  in  highly  forward  scattering  random  media. 

Tasks  Completed: 

The  tasks  completed  during  this  period  follow. 

Scalar  wave  diffusion  equation(SWDE),  was  derived  to  describe  the  non-diffusive 
effects  taking  place  at  large  distances  from  a  source.  SWDE  has  been  shown  to  be 
superior  to  other  available  modifications  of  the  diffusion  equation  to  describe  non- 
diffusive  effects  in  isotropically  scattering  media  (publication  [3]). 

The  new  concept  of  photon  paths  in  highly  forward  scattering  media  was 
investigated  using  time-resolved  experiments  and  the  theory  derived  from  non-Euclidean 
diffusion  equation  (NED).  NED-based  concepts  of  average  and  most  probable  (“Fermat”) 
photon  paths  were  successfully  used  to  interpret  the  experimental  results  which 
demonstrate  in  contrast  to  the  predictions  of  the  diffusion  theory,  that  photons  “prefer” 
certain  paths  on  a  scale  on  the  order  of  transport  mean  free  paths  (publications  [1],  [4]). 

Time-and-angle-resolved  experiments  have  demonstrated  that  the  non-diffusive 
effects  in  photon  migration  are  appreciable  even  at  distances  as  large  as  10-151,  (20-30 
mm)  from  the  source  (publications  [2],  [5],  [6]). 
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Comparison  of  the  above  time  resolved  data  with  the  results  based  on  NED  has 
demonstrated  the  potential  of  NED  to  describe  the  nondiffusive  multiple  scattering  effect 
on  photon  migration  in  turbid  media  (publications  [2],  [6],  [7]). 

Results: 


The  following  highlights  our  results: 

Photon  density  modes  beyond  the  diffusion  approximation  in  a  turbid  medium 
with  isotropic  scatterers. 


Various  modifications  of  the  diffusion  model,  available  to  date,  designed  to 
describe  photon  migration  in  scattering  and  absorbing  media  are  summarized  in  the 
following  equation: 


a- 


d2N 


.  dN 


dt: 


— +  (v  +  pvj— A/V+va(v+yvJV  =  0 


(1) 


for  the  photon  number  density  N(r,  t) .  Here  v  and  vfl  are  the  transport  and  absorption 

collision  frequencies.  Different  sets  of  numerical  parameters  {a,p,y}  corresponding  to 
different  models  can  be  found  in  the  literature:  {0,0,0}  (Furutsu,  1994),  {0,1,1} 
(Paterson,  Chance,  Wilson,  1989),  {1,2,1}  (Ishimaru,  1989,  Kashke,  1993)  and  {1/3, 
2/3, 1/3}  (Lax,  1993,  Morse  &Feshbach,  1963).  Starting  with  the  general  radiative  transfer 
equation  we  have  shown  than  none  of  the  above  versions  of  Eq.(l)  is  actually  valid  when 
used  in  non-diffusive  domain,  and  a  better  new  equation,  called  Scalar  wave  Diffusion 
Equation  (SWDE)  was  derived.  SWDE  coincides  in  form  with  Eq.(l),  but  the  set  of 
coefficients  {a,p,y}  is  totally  different:  {a,p,y}  =  { l/5,2/5,l/5}.  To  check  the  available 
versions  of  Eq.(l)  and  SWDE,  analytical  solutions  of  these  equations  N  were  compared 
with  a  numerical  solution  NE  of  the  exact  radiative  transfer  equation  for  the  intensity 
modulated  source  T|  switched  on  at  a  time  t=0  : 

r|  =  ~Q(t)S(r)exp(  im)  (2) 

47C 


Relative  changes  in  amplitude  *P  =  (jw|-|/V£|)/|V£|  are  plotted  versus  co/v  for  source 

detector  distance  r  =  25  1,  (Fig  la)  and  r  =  50  1,  (Fig  lb).  The  parameter  \j t  shows  how 
close  an  approximate  solution  is  to  the  reference  NE  .  As  it  is  clear  from  Fig  1,  the  SWDE 
solution  is  much  closer  to  the  exact  numerical  solution,  NE,  than  any  other  model. 

Note  that  the  absorption  collision  frequency  va  is  completely  equivalent  to  the 
parameter  p  in  the  Laplace  transform  of  the  experimental  temporal  profiles: 

oo 

I(p)  =  j e~pt I (t)dt .  When  p  =  id)  is  an  imaginary  number,  this  corresponds  to  an 
o 
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intensity  modulated  probing  photon  beam.  High  values  of  |p|=0.1v  select  early  arriving 

light  and  correspond  to  smaller  blurring.  Therefore  SWDE  may  lead  to  an  improved 
resolution  when  used  in  early-light-based  tomography. 

A  generalization  of  SWDE  to  the  case  of  anisotropic  scattering  modes  is  to  be 
considered  separately. 

Theory  and  experiment  on  ultrashort  pulse  non-diffusive  propagation  in  highly  forward 
scattering  media. 

Time-and-angle  resolved  intensity  measurements  of  scattered  light  were  using 
femtosecond  laser  pulses  at  625  nm  in  a  highly  forward  scattering  medium  at  various 
source-detector  distances  much  larger  than  1(.  The  detector  was  placed  in  the  direction  of 
initial  laser  beam  propagation,  and  orientation  of  the  fiber  was  altered  (Fig  2a,  inset).  The 
medium  was  0.16%  (by  weight)  suspension  of  1.01  |im  polystyrene  spheres  in  water  with 
1,  =2 mm  ,  la  =  300mm,  g  =  0.93. 

Fig.  2  shows  the  experimental  results  for  various  source  detector  distances.  Note  the 
absence  of  ballistic  component  for  highly  forward  scattering  media.  Early  arriving  light 
is  represented  by  snake  photons.  One  may  note  that  even  at  distances  as  large  as  10-15  lt 
(20-30  mm)  angular  distribution  of  light  is  much  too  far  from  being  isotropic.  Highly 
isotropic  light  distribution  is  a  condition  required  for  the  diffusion  approximation  to  be 
valid.  These  results  show  that  diffusion  approximation  (DA)  is  hardly  a  reliable  tool  for 
high  resolution  optical  tomography.  The  above  experimental  data  were  used  to  check  on 
predictions  based  on  non-Euclidean  diffusion  equation  (NED),  a  possible  alternate  to  DA. 

Approximate  solution  to  the  non-Euclidean  diffusion  equation 

During  the  reported  period  we  have  found  approximate  analytic  solution  to  the  NED: 

+  sVrn  -  DsAsn + v an  =  8(tj8fr)5(s  -  s0 )  (3) 

at 

for  the  specific  photon  number  density  n(t,r,s)  as  a  function  of  time  t,  position  r,  and 
velocity  of  photons  s.  In  NED,  Eq  (3)  scattering  is  described  in  terms  of  photon  diffusion 
on  the  velocity  sphere,  where  Ds=c4/6D  is  the  photon  diffusion  coefficient  on  the  velocity 
sphere  expressed  in  terms  of  the  ordinary  3D-diffusion  coefficient  D.  We  have  found 
exact  path  integral  solution  of  Eq(3)  and  analytically  calculated  this  path  integral  to 
obtain  the  photon  number  density  N(r,  t)  =  J  dsn( t,r,s)  using  cumulant  expansion  in 
random  process  theory: 

A Hr,  0  =  7  1  . =  * expl—  ~  A~  p  (r  -  rc  )a  (r  -  rc)p  L_Va'  (4) 

iJ(4k)3  det  Aap  l  4  PJ 

In  Eq.(4),  rc(t)  =  s0(ll/c)fl 
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\#<soO  =  £*#fy-fi  +^f3>  +  -£i<’M*o>p<f<  ~\h -A2). 

ct 

X=— ,  fn  =\-exp(-n%). 

h 


Eq  4  describes  ultrashort  pulse  propagation  represented  as  translational  motion  and 
spread  of  an  initially  point-like  photon  cloud  launched  in  s0  direction.  In  Fig.3  the 

solution  N(r,t)  to  Eq  (4)  is  compared  with  the  experimental  results  of  Fig.2.  The 
corresponding  angular-dependent  experimental  profiles  were  integrated  over  all  detector 
orientations  to  result  in  experimentally  estimated  photon  number  density.  The 
experimental  results  confirm  superiority  of  NED  to  DA  in  describing  the  early  arriving 
light  intensity  and  the  whole  temporal  profile.  Note  that  the  approximate  solution  takes 
causality  principle  into  account  qualitatively.  Obtaining  analytical  results  for  NED-based 
specific  number  density  n(  t,r,s)  and  its  direct  comparison  with  experimental  data  are 
under  way.  We  also  plan  to  find  exact  analytic  representation  for  early  portion  of 
scattered  light  intensity  with  the  causality  principle  taken  into  account  exactly.  These 
results  may  be  used  as  a  basis  for  future  early-arriving-light  tomography. 


NED  and  nondiffusive  structure  of  photon  path  in  random  media. 

To  examine  the  photon  path  pattern,  time  resolved  intensity  measurements  of  scattered 
light  for  a  fixed  source(S)-detector(D)  arrangement  have  been  used  to  determine  the 
position  at  which  a  small  absorbing  object  equidistant  from  the  source  and  detector 
produces  maximum  effect  on  the  detectors  readings.  Launching  and  collecting  fibers  were 
parallel  to  each  other  and  perpendicular  to  SD  line.  The  experiment  revealed  in  contrast  to 
the  predictions  of  the  diffusion  theory,  that  photons  traverse  certain  curvilinear  paths 
between  the  source  and  the  detector.  These  results  can  be  interpreted  in  terms  of  average 
photon  trajectory  ra(t\T)  between  a  source  and  a  detector,  corresponding  to  a  certain 

total  travel  time  T.  Note  that  ra(t\T)  in  diffusion  approximation  is  an  exactly  straight 
line  between  a  source  and  detector.  The  point  on  ra(t\  T)  for  t  =  T/2  is  the  equdistant 
from  the  source  and  the  detector. 

We  have  used  the  analytical  solution  for  the  NED,  Eq(4),  to  estimate  ra(t\T): 

(ra(t\T))i  =  (A(s0,t)  +  A(sf0,T-t))~/  *(Ajk(sf  T-t$( t)sok  +  Ajlc (s0 t)( -Q(T-t)sJk)) 

(5) 

Here  r  is  the  position  of  the  detector  with  respect  to  the  source,.  J0  ,Sf  are  the  initial 
and  final  velocities  of  photons  (sj=  -s0),  Q(t)  -  (l,  /c)f{(t),  and  functions  //tj 
and  A^fjo/jwere  determined  in  connection  with  Eq.(4). 
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Fig  4.  shows  calculated  average  photon  paths  corresponding  to  the  source-detector 
distance  R=22mm  and  travel  time  T=157ps  for  various  values  of  1, .  The  value  1,  =  15mm 
corresponds  to  the  condition  of  our  experiment.  The  cross  (x)  shows  the  experimentally 
determined  position  of  the  object  corresponding  to  the  maximal  decrement  in  the  detector 
reading.  Satisfactory  agreement  of  the  theory  and  the  experiment  can  be  noted.  When  the 
ratio  1/L  is  not  small  ra(t  I T)  coincide  with  the  most  probable  (Fermat  paths  rF  ).  As  the 

ratio  1/T  decreases  the  population  of  path  increases  and  ra  deviate  from  rF  toward  SD  line 
in  accordance  with  experimental  observation  of  smaller  lt. 
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ACCOMPLISHMENTS: 

1.  Introduced  and  tested  was  a  new  scalar-wave-diffusion  equation(SWDE)  aimed  at 
describing  nondiffusion  effects  in  a  turbid  media  with  isotropic  scatterers  for  frequency 
modulation  and  time  resolved  studies. 

2.  Time  and  angle  resolved  experiments  demonstrated  that  non-diffusive  effects  in 
photon  migration  are  appreciable  at  distances  as  large  as  10-15  1,  (20-30  mm).  These 
effects  promise  to  be  important  for  early  light  tomography  applications. 

3.  Comparison  of  time-and-angle-resolved  experiment  with  the  non-Euclidean  diffusion 
equation(NED)  introduced  earlier  showed  that  NED  can  be  regarded  as  a  feasible 
approach  to  describe  nondiffusive  photon  path  patterns  in  highly  forward  scattering 
media. 
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Advances  in  Optical  Imaging  and  Photon  Migration,  2,  14-17 
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Interaction  with  optical  imaging  community 

Quotation  of  our  papers  published  in  the  frame  of  this  grant  in  1996: 

1.  L.Yang,  J.U.Beecham  1996:  Time  resolved  Anisotrophy  Spectroscopy  and 
Imaging  theory  Scattering  Media  using  Most  Favorable  Fermat  Photon  Paths, Biophysical 
Journal,338.  (the  idea  of  most  probable  photon  paths  was  used  to  interpret  experimental 
results) 

2.  S.B.Colak,  H.Schomberg,  G.W.’t  Hooft,  M.B.Van  der  Mark,  1996:  Optical  Back 
Projection  Tomography  is  Heterogeneous  Diffusive  media,  OSA  TOPS  on  advances  in 
optical  Imaging  and  Photon  Migration,  V2,  282.  (the  idea  of  most  probable  photon  paths 
was  discussed) 


The  following  groups  requested  information  about  our  work  done  under  this 

grant: 


1. 

2. 

3. 

4. 

5. 

6. 


George  Weiss 
Mark  Brezinski 
Michael  Feld 
Norman  McMormick 
Arjun  Yodh 
Steve  Jacques 


( National  Institute  of  Health) 

( Massachusetts  general  Hospital) 

( Massachusetts  Institute  of  Technology) 
( University  of  Washington) 

( University  of  Pennsylvania) 

( Oregon  Medical  Laser  Center) 
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Figure  Captions 


Fig.  1  Comparison  of  the  numerical  solution  of  the  radiative  transfer  equation  NE  with  the 
solution  N,  corresponding  to  the  scalar-wave-diffusion  equation  and  the  different 

modification  of  Eq(l).  Relative  change  in  amplitude  'F  =  f|A^|-|A^£|j/|A^£|  is  plotted:  0  - 

exact  numerical  solution  NE,  1-SWDE,  2-Ishimaru,  1989,  Kashke,  1993,  3-Paterson, 
Chance,  Wilson,  1989,  Fishkin,  Gratton,  1993;  4-Furutsu,  1994;  5-Lax,  1993,  Morse  & 
Feshbach,  1963.  2’  corresponds  to  the  curve  2  scaled  as  1:2500.  (a)  source  detector 
distance  r  =  251t,  va  =  0;  (b)  r  =501t,  va  =  0. 1  v. 

Fig2.  Intensity  profiles  measured  at  various  source  detector  distances  r  and  for  various 
orientations  of  the  detecting  fiber:  (a)  r  =  51t,  (b)  r  =  71t,  (c)  r  =  10  It,  and  (d)  r  =  151t,  The 
inset  shows  the  geometry  of  the  experiment. 

Fig.3  (a)  Comparison  of  the  measured  and  calculated  NED-based  and  DA-based  photon 
number  density  profiles  N(t)  at  n=51,;  (b)  Full  width  half  maximal  of  N(t)  and  various 
distances;  (c)  Tp  position  of  the  maximum  of  the  N(t)  curve;  (d)  Ip(r)/I(r=151t)  relative 
amplitude  of  the  N(t)  curve.  The  most  probable  (Fermat)  path  is  indicated  by  symbol  F. 

Fig  4  Calculated  average  photon  paths  corresponding  to  the  source-detector  distance 
R=22  mm  and  travel  time  T=157ps  (  overall  path  length  35  mm)  for  various  values  of  1, : 
(1)  15  mm;  (2)  9  mm;  (3)  4.5  mm;  (4)  0.5  mm.  The  cross  corresponds  to  the  experimental 
data  taken  at  1=  15mm. 
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STATISTICS: 

1 )  Number  of  papers  published  in  refereed  journals  3 

2)  Number  of  papers  submitted  or  in  press,  refereed  journals  1 

3)  Number  of  books  or  chapters  published,  refereed  non-serial  publications  0 

4)  Number  of  books  or  chapters  submitted  or  in  press,  refereed  non-serial  0 

5)  Number  of  invited  presentations  at  scientific  conferences  1 

6)  Number  of  contributed  presentations  at  scientific  conferences  3 

7)  Number  of  technical  reports  and  papers  in  non-refereed  journals  2 

8)  Number  of  undergraduate  students  supported  (at  least  part  time)*  0 

9)  Number  of  graduate  students  supported  (at  least  part  time)*  0 

10)  Number  of  post-docs  supported  (at  least  part  time)*  2 

1 1)  Number  of  other  professional  personnel  supported  (at  least  part  time)*  0 

12)  Number  of  female  grad  students  0 

13)  Number  of  minority  grad  students  2 

14)  Number  of  Asian  grad  students  1 

15)  Number  of  female  post-docs  0 

16)  Number  of  minority  post-docs  1 

17)  Number  of  Asian  post-docs  1 

*  Part  time’  means  that  at  least  25%  support  this  year  on  this  ONR  grant.  NOTE 
Minorities  include  Blacks,  Aleuts,  American  Indians  and  Hispanics  only. 
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